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Abstract. This paper develops strong solutions and stochastic solutions for the 
tempered fractional diffusion equation on bounded domains. First the eigenvalue 
problem for tempered fractional derivatives is solved. Then a separation of variables, 
and eigenfunction expansions in time and space, are used to write strong solutions. 
Finally, stochastic solutions arc written in terms of an inverse subordinator. 



1. Introduction 

Transient anomalous sub- diffusion equations replace the first time derivative by a 
tempered fractional derivative of order < /3 < 1 to model delays between movements 
[H [7] . These governing equations have proven useful in finance [21 [6] and geophysics 
[TTj to model phenomena that eventually transition to Gaussian behavior. A sto- 
chastic model for transient anomalous sub-diffusion replaces the time variable in a 
diffusion by an independent inverse tempered stable subordinator. This time-changed 
process is useful for particle tracking, a superior numerical method in the presence 
of irregular boundaries [221 123]. The idea of tempering was introduced by Mantegna 
and Stanley [131 HI] and developed further by Rosinski [T9] . 

Section [2] provides some background on diffusion and fractional calculus, to estab- 
lish notation, and to make the paper relatively self-contained. Section [3] uses Laplace 
transforms and complex analysis to prove strong solutions to the eigenvalue problem 
for the tempered fractional derivative operator. Then Section H] solves the tempered 
fractional diffusion equation on bounded domains. Separation of variables and eigen- 
value expansions in space and time lead to explicit strong solutions in series form. 
Stochastic solutions are then developed, using an inverse tempered stable time change 
in the underlying diffusion process. 
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2. Traditional and fractional diffusion 

Suppose that D is a bounded domain in M'^. A uniformly elliptic operator in 
divergence form is defined for u G C'^{D) by 

(2.1) Ldu= y]--{a,,{x)p-] 



with aij{x) = aji{x) such that for some A > we have 



nd 



(2.2) A ^ < J2 ^iA^)yiyo < A"' ' y ^ 

We also assume that for some A > we have 

n 

(2.3) Y\aij{x)\<k forallxeD. 

Take a = aa'^ , and B{t) a Brownian motion. Let X{t) solve the stochastic differential 
equation dX{t) = b{X{t))dt + a{X(t))dB(t), and define the first exit time td{X) = 
mi{t > : X{t) ^ D}. An application of the Ito formula shows that the semigroup 

(2.4) TD{t)f{x) = E,[f{X{t))I{TD{X) > t)] 

has generator (12. ip . see Bass [H Chapters 1 and 5]. Since TD{t) is intrinsically ul- 
tracontractive (see [H Corollary 3.2.8, Theorems 2.1.4, 2.3.6, 4.2.4 and Note 4.6.10] 
and [lOl Theorems 8.37 and 8.38]), there exist eigenvalues < ?7i < 772 < 773 ■ ■ ■ , with 
rjn ^ 00 and a complete orthonormal basis of eigenfunctions ipn in L^{D) satisfying 

(2.5) Loipnix) = -r]ni^n{x), X E D : tpnldD = 0. 

Then p£){t,x,y) = J2'^=i^~^"^'^n{x)ipn{y) is the heat kernel of the killed semigroup 
Tp. This series converges absolutely and uniformly on [tg, 00) x D x D for all to > 0. 
Denote the Laplace transform (LT) t — )■ s of u{t,x) by 

POO 

u{s, x) = Ct[u{t, x)] = / e~^^u{t,x)dt. 

Jo 

The 'j/'n-transform is defined by u(t, n) = il)n{x)u{t, x)dx and the ^/'n-Laplace trans- 
form is defined by 

(2.6) u{s,n) = I il)n{x)u{s,x)dx. 

Jd 

Since {ipn} is a complete orthonormal basis for L^{D), we can invert the -j/^^-transform 
to obtain u(t,x) = '^^ u(t,n)ipn{x) for any t > 0, where the series converges in the 
L2 sense [201 Proposition 10.8.27]. 
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Suppose that D satisfies a uniform exterior cone condition, so that then each x e 
dD is regular for [3], Proposition 1, p. 89]. If / is continuous on Z), then 

u{t,x) = Tn{t)f{x) = E,[f{X{t))I{Tn{X) > t)] 
^^■^^ = / PD{t,x,y)f{y)dy = f2e-^-'M^)f{n) 

solves the Dirichlet initial-boundary value problem [HI Theorem 2.1.4]: 

du{t, x) 



, , dt 
(2.8) 



L[)u{t,x), X G -D, t > 0, 



u{t, x) = 0, X e dD, 
u{0,x) = /(x), X E D. 

This shows that the diffusion X{t) killed at the boundary dD is the stochastic solution 
to the diffusion equation (12. 8p on the bounded domain D. 

The Riemann-Liouville fractional derivative [T8l 121] is defined by 

(2.9) tMt)- ' ^r'^M^' 



dt^"' ' r(i-/3)rft io {^-^y 

when < /3 < 1. The Caputo fractional derivative [5] is defined by 

(2.10) ' 



.ay r(i-/3)i„ (t-»)« 

when < /3 < 1. It is easy to check using £[t~^] = s^~^V{\ — f3) that 



(2.11) £ 



while the Caputo fractional derivative (I2.1QP has Laplace transform s^g{s) — s^ ^di^)- 
It follows that 



dt^"^ ' \dtj ' T{l-P)' 

Substituting a Caputo fractional derivative of order < /3 < 1 for the first order 
time derivative in (12. 8 p yields a fractional Cauchy problem. This fractional diffusion 
equation was solved in [151 Theorem 3.6]. Those solutions exhibit anomalous sub- 
diffusion at all times, with a plume spreading rate that is significantly slower than 
(12. 7p . Many practical problems exhibit transient sub-diffusion, converging to (12. 7p at 
late time [21 |6l [17]. Hence, the goal of this paper is to extend the results of [15] to 
transient sub-diffusions. 
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3. Eigenvalues for tempered fractional derivatives 



Suppose D{x) is a standard stable subordinator with Levy measure (j){y, oo) = 
I{y > 0)?/-^/r(l - /3) for < /3 < 1 , so that E[e-"-°(^)] = e'^^^^), where the 
Laplace symbol ^{s) = = /o°° (1 ~ ^~^^) (l>i3{dy). If fx{t) is the density of D{x), 



then q\(t,x) = fx{t)e ^^/e is a density on x > with Laplace transform (LT) 



(3.1) 



qx{t,x) dt 



Ut)dt 



-xipx{s) 



where ip\{s) = (s+A)^ — A''. Rosinski [19] notes that the tempered stable subordinator 
D\{x) with this Laplace symbol has Levy measure (l)\{dy) = e'^^ (p (i{dy) . 

Define the Riemann-Liouville tempered fractional derivative of order < /3 < 1 by 



(3.2) 



■git) 



1 d fe^'g{s)ds 



r(i - p) dt Jo [t - s)/3 

as in [1]. We say that a function is a mild solution to a pseudo-differential equation 
if its Laplace transform solves the corresponding equation in transform space. 



Proposition 3.1. The density q\{t,x) of the tempered stable subordinator D\{x) is 
a mild solution to 

qx{t,x). 



(3.3) 

Proof. Clearly (-5, x) 
(3.4) 



^gA(t,x) 



e solves 
d 



dx 



q\{s,x) = -ipxis) qx{s,x) 



with initial condition qx{s,0) = 1. The right-hand side of ( 13. 4 p involves a pseudo- 
differential operator ipx{dt) with Laplace symbol ipx{s), see Jacob [11]. To complete 
the proof, it suffices to show that ipx{s)g{s) is the LT of (13. 2p . Since C[e'^^g{t)] = 
g{s - A), we get 

d^ 



(3.5) 
which leads to 
(3.6) 



C 



dtp 



{e''git)) 



s^gis 



e-"^(e^(t)) 



A), 



Then (13. 3p follows easily. This also shows that ipx{dt) is the negative generator of the 
Co semigroup associated with the tempered stable process. □ 

Define the inverse tempered stable subordinator 



(3.7) 



Exit) = mf{x > : Dxix) > t}. 
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A general result on hitting times [111 Theorem 3.1] shows that, for all t > 0, the 
random variable Ex{t) has Lebesgue density 

(3.8) g^{t,x)=l (f)x{t - y,oo)qx{y,x)dy 

Jo 

and (t,x) I— gx(t,x) is measurable. Following [121 Remark 4.8], we define the Caputo 
tempered fractional derivative of order < /3 < 1 by 

Proposition 3.2. The density fl3.8p of the inverse tempered stable suhordinator (13.71) 
is a mild solution to 

d \ ^'"^ 

(3-10) di^'^^'''^ = ~[m) 

Proof. Theorem 4.1 in [16] shows that (13. 8p is a mild solution to the pseudo-differential 
equation 

d 

(3.11) — gx{t, x) = -i)x{dt)g\{t, x) + 5(x)0A(t, oo) 

where ipx{dt) is the pseudo-differential operator with Laplace symbol ipx{s), i.e., the 
Riemann-Liouville tempered fractional derivative (13. 2p . Use 

1 /"°° 

(3.12) 0,(t,oo) = ^;^^-^y^ e-'^(3r-P-'dr 
to rewrite (13. lip in the form 

(3.13) = _^,,(,,) + 

and then apply I^M) with 5fA(0,x) = 6{x) to get fIXTU]) . □ 

The next two results establish eigenvalues for Caputo tempered fractional deriva- 
tives, which will then be used in Section H] to solve tempered fractional diffusion 
equations by an eigenvalue expansion. 

Lemma 3.3. For any fi > 0, the Laplace transform 

POO 

(3.14) gxit,fi) = C4gxit,x)]= / e-'^^gxit^x) dx = E[e-''''^^'^] 

Jo 

is a mild solution to the eigenvalue problem 

(3.15) — g\{t,fi) = -fig\{t,fi) 



dtj 

with gx{0,fi) = 1, for the Caputo tempered fractional derivative (13.90 . 



Proof. Equation (3.12) in [TB] shows that 

(3.16) Ct[Mt,oo)] = s-'Ms). 

Then (13. 8p together with the LT convolution property shows that 

(3.17) = Ct[gx{t,x)] = Ct[Mt,oo)]Ct[qx{t,x)] = -^^p^{s)e-^^'^'^ 

for any x > 0, and then a Fubini argument shows that the double Laplace transform 



Ms) 



-(/i+i/'A(s))3::^^ 



(3.18) Gx{s,fi) = CtC4gx{t,x)] 
Rearrange (13.181) to get 

(3.19) - fiGxi-s, /i) = Ms)Gxis, /i) - s-Va(s), 
use (13.161) along with (13. 9p and (13.121) to see that 



Ms) 



s{^ + ^Ijx{s))' 



(3.20) 



^\{s)Gx{s,n) - s Va(s), 



then substitute into (I3.19P to get 



This proves that (13.141) is the mild solution to (I3.15p . 



□ 



The next theorem is the main technical result of this paper. It shows that Laplace 
transforms (I3.14p of inverse tempered stable densities are the eigenvalues of the Ca- 
puto tempered fractional derivative, in the strong sense. 



Theorem 3.4. ForO < (3 <1, let 



-tx 



(3.21) 



k{t) 



vr sin(/37r) 

For any yU, A > 0, 7^ A^, the function gx{t, jj) in (I3.14p can he written in the form 

POO 

(3.22) gxit,^I) = ^ / (r + A)-^e-*(^+^)$(r, l)dr, 



where 



$(r,r 



sin(/37r) 



r2/3 sin2(/37r) + {fi - + cos(/37r))2 
and 

(3.23) \dtgx{t,fi)\<fik{t). 

Then gx(t,fi) is a strong (classical) solution of the eigenvalue problem (I3.15p . 



Proof. The proof extends Theorem 2.3 in Kochubei [12] using some probabihstic argu- 
ments. Since Exit) has continuous sample paths, a dominated convergence argument 
shows that gx{t, jj) = E[e~'^-^^'^*^] is a continuous function of t > 0. Use fl3.18p to write 

(3.24) G.(.. ^) = )1- ''^W - + 



and note that Gx{s,n) is analytic off the branch cut arg(s) = tt, \s\ > 0. The Laplace 
inversion formula [9l p. 25] shows that for suitable 7 > and for almost all t > 0, 

3.25 g^[t,fi) = —— — — — -ds. 

at 2txi J^^ioo s [{s + X)p - Ap + /ij 

Let \ < oj < 1 and consider the closed curve C-.^^^^ in C, formed by a circle of radius 
with a counterclockwise orientation, cut off on the right side by the line Im [z) = 7, 
and by the curve S^^^j on the left side, consisting of the arc 



and the two rays 



T-y^uj = {s E C : |s| = 7, I args| < un} 

r^,^^ = e C : |s| > 7, |args|=W7r}, 
^y,oj = {-s^C- |s|>7, I args| = -wtt}. 
By Cauchy's Theorem, the integral 



ds = 



and then Jordan's Lemma P, p. 27] implies that we can let i?„ — )■ 00 to get 

3.26 g\(tn) = - / —-^ ^L^ds. 

Now pass the derivative inside the integral to get 

which also implies the smoothness of the function t — )■ gx(t,fi). It is not hard to 
check that the integral over T^^^ tends to zero as 7 ^ 0. Compute the remaining 
path integral, and let a; — 1 to obtain (13. 22 p . 
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Differentiate (I3.22p with respect to t and use r'' sin(/37r)$(r, 1) < 1 to write 



\dtgx{t,n)\ 



(3.28) 



< 



TT sin(/37r 



r + A)-i[afe-*("+^)]$(r,l)dr 
t^^ir(l -/3) = /iA;(t), 



7rsin(/37r 

so that (13.231) holds. Note that \gx(t,fi)\ < 1, and write 
d 



< 



1 

r(i-/3) Jo 
1 



r(i - /3) 7o 
r(i-/3)/ 



As ^[£'A(s,/i)] \ 



Ae'''^|£rA(s,/i)| + e 

As , 



Ae 



7rsin(/37r) 



9s 



(t - s)/^ 
(is 



it 



Then a simple dominated convergence argument shows that the Riemann-Liouville 
fractional derivative of e^^g\{t,^) is a continuous function of t > 0. Now it follows 
from (13. 2p and (13. 9p that the Caputo tempered fractional derivative of gx{t,fi) is 
continuous in t > 0. Since both sides of (I3.15P are continuous in t > 0, it follows from 
Lemma 13.31 and the uniqueness theorem for the Laplace transform that (13.151) holds 
pointwise in t > for all /i > 0. □ 



4. Tempered fractional diffusion 

Let = (0,oo) X D, define Hl^{D^) = {u: D^^^: LDu{t,x) E 0(0^)}, 
and let nijD^) = Hl^{D^) n {u : \dtu{t,x)\ < k{t)g{x), g G L^{D), t > 0}, 
where k{t) is defined in (13.211) . 

Theorem 4.1. Let D be a bounded domain with dD G C^'" for some < a < 1, and 
let X(t) be a continuous Markov process with generator (12.10 . where aij G C°'{D). 
Then, for any f G D{L£,) H C^{D) fl C'^{D) such that the eigenfunction expansion 
of Lof with respect to the complete orthonormal basis {ipn} converges uniformly and 
absolutely, the (classical) solution of 

L£)u(t,x), X E D, t > 0; 
0, X e dD, t> 0; 

/(x), xeD, 
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(4.1) (^) uit,x) ^ 

u(t,x) = 
u{0,x) = 



for u e ni^{D^) n Cb{D^) n C\D), is given by 

u{t,x) = E,[f{X{E^mHrD{X) > = E,[/(X(E,(t)))/(r,,(X(E,)) > t)] 

poo 

(4.2) = / TD{l)fix)gx{tJ)dl = y2fi^)M^)9xit,Vn), 
Jo 

where Ex{t) is defined by (13. 7p . independent of X{t), gx{t,r]) = E(e~''-^^(*^) is its 
Laplace transform, and T^it) is the killed semigroup (12.41) . 

Proof. The proof uses a separation of variables. Suppose u(t,x) = G(t)F{x) is a 
solution of (14.11) . substitute into (14. ip to get 



d 



F{x) [j^j G{t) = G{t)LnF{x) 
and divide both sides by G{t)F{x) to obtain 

iif'Git) _LnF{x) 



Git) F{x) 



-T]. 



Then we have 



d 



/3,A 



(4.3) (^-J G{t) = -r^Git), t > 
and 

(4.4) LF{x) = -r]F{x), x e D, F\dD = 0. 

The eigenvalue problem (14.40 is solved by an infinite sequence of pairs {(/i„, ■?/;„)}, 
where < ?7i < 772 < 773 < ■ ■ ■ , ?7„ — )■ oo, as n — )■ oo, and ipn forms a complete 
orthonormal set in L^{D). In particular, the initial function / regarded as an element 
of L'^{D) can be represented as 



(4.5) /(a;) = ^/»V;, 



, X) 
n=l 



Use Lemma [373] to see that Gn{t) = f{n)gx{t,rin) solves (14. 3p . Sum these solutions 

i;n{x)Gn{t) to dHD, to get 

00 

(4.6) u{t, x) = ^ f{n)gx{t, Vn)^n{x). 

n=l 

It remains to show that (14.61) solves (14. ip and satisfies the conditions of Theorem 14.11 
The remainder of the proof is similar to [151 Theorem 3.1], so we only sketch the 
argument. First note that (14.60 converges uniformly in t G [0, 00) in the sense. 
Next argue ||M(t, ■) — f\\2,D — as t — using the fact that, since gx(t,X) is the 
Laplace transform of Ex{t), it is completely monotone and non-increasing in A > 0. 
Use the Parseval identity, the fact that is increasing in n, and the fact that gx{t, fJ^n) 
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is non-increasing in n > 1, to get ||u(t, ■)||2,d < ^7a(^5/^i)||/||2,_d- A Fubini argument, 
which can be rigorously justified using the bound \dtu(t,x)\ < k{t)g{x) from Theorem 
I3.4[ shows that the ipn transform commutes with the Caputo tempered fractional 
derivative. For this, it suffices to show that the '?/'„-transform commutes with the 
Caputo fractional derivative of e^^u(t,x). To check this, write 



D 



at) 



{e^^u{t,x)) dx 



D 



r(i-/3) 



9 (e'^*u(s, x)) ds 



ds 



it - sr 



dx 



r(i - /3) Jo 

1 



d 



D 



'^ni.x)— (e^'*M(s,x)) dx 



ds 



r(i - 13) Jo ds 



* d 

— (e^^M(s,n)) 



tljn{x)u{s, x)dx 

ds 



it - sr 

ds 



(t - sY 



it-s) 







Then the Caputo tempered fractional time derivative and the generator Ld can be 
applied term by term in (14.61) . Next show that the series (14. 6 p is the classical solution 
to (14. ip by checking uniform and absolute convergence. Argue that u G C^{D) using 
p!Ol Theorem 8.33], and the absolute and uniform convergence of the series defining /. 
Finally, obtain the stochastic solution by inverting the T/Cn-Laplace transform. Since 
{il^n} forms a complete orthonormal basis for L'^[D), the '^^-transform of the killed 
semigroup TD{t)f{x) = J2m=i ^'^""^i^mix) f (m) from (12.70 is given by 



(4.7) 



[Tn{t)f]{n) = e-"^-f{n). 



Use Fubini together with ( 14. 6p and (14. 7p to get 



u{t,x) = ^/(n)V^„(x)^A(t,/in) = J^V'nla:) / f\n)e ^-ygx{t,y)dy 

= y'^/'nla;) / [TDiy)f]{n)gx{t,y)dy 
n=i -^0 



n=l 



gx{t,y)dy 



TD{y)f{x)g\{t,y)dy 



E,[/(X(E,(t)))/(r^(X) > Exit))]- 
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The argument that 



E,[f{X{E^mHrD{X) > E,m = E,[f{X{E^mHrD{X{E^)) > t)] 

is similar to [151 Corollary 3.2]. Uniqueness follows by considering two solutions ui, U2 
with the same initial data, and showing that Ui — ^2 = 0. □ 

Remark 4.2. In the special case where = A, the Laplacian operator, sufficient 
conditions for existence of strong solutions to (14.1 p can be obtained from [T^ Corollary 
3.4]. Let / G Cf^i^D) be a 2/c-times continuously different iable function of compact 
support in D. If /c > l+3(i/4, then (14. ip has a classical (strong) solution. In particular, 
if / e C^{D), then the solution of dH]) is in C°°{D). 

Remark 4.3. In the special case where L/) = A on an interval (0, M) C M, eigenfunc- 
tions and eigenvalues are explicitly known, and solutions to the tempered fractional 
Cauchy problem can be made explicit. Eigenvalues of the Laplacian on (0, M) are 
(mr/M)'^ for 77, = 1, 2, ■ ■ ■ and the corresponding eigenfunctions are -^sm{mrx/M). 
Using this eigenfunction expansion, the solution reads 

oo oo 

u{t, x) = ^ f{n)ipnix)gx{t, = ^ f{n) sm{mTx/M)gx{t, {rni/Mf). 

n=l n=l 
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